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ABSTRACT 

It is pointed out that the beaming correction commonly inferred from the 
achromatic breaks in the afterglow light curve, is biased in situations where 
the isotropic equivalent energy is affected by factors other than the spread in 
opening angles. In particular, it underestimates the beaming factor of sources 
observed off-axis. Here we show that both the slopes and scatters in the Am- 
ati and Ghirlanda relations, and the difference between them, are quantitatively 
consistent with a model proposed recently, in which the E iso — v peak relation, as 
originally derived by Amati et al., is due to viewing angle effects. The quantita- 
tive difference between them confirms the relations between opening angle and 
break time suggested by Frail et al. 

Subject headings: black hole physics — gamma-rays: bursts and theory 

A correlation between the isotropic equivalent energy of GRB's, E iso , and the location of 
the spectral peak at local redshift, z/ p , has been reported recently for a sample of BeppoSAX 
and HETE2 sources with measured redshifts and well observed spectra (Amati et al.2002; 
Atteia et al. 2003; Lamb et al. 2004). This relation can be represented as a power law, 

(E iso /10 52 erg) ~ (hz/p/lOOKeVr, (1) 

with a = 2, and spans two orders of magnitude in u p , roughly from 10 KeV to 1 MeV, and 
about four orders of magnitude in Ei so . The extension to low energies is based on a small 
number of X-ray flashes with measurable redshifts. It does not preclude the possibility that 
many X-ray flashes derive their low apparent peak luminosities in part due to cosmological 
redshift. A similar correlation between the collimation-corrected energy and u p , with a 
smaller scatter and a different slope, has been reported by Ghirlanda et al. (2004; hereafter 
GGL04) for a subsample of the sources that exhibit achromatic breaks in the afterglow light 
curves. 
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Eichler & Levinson (2004; hereafter EL04) proposed that the E iso — v v correlation can be 
interpreted as due to an orientation effect. They have shown that an annular jet with a single, 
universal spectrum can give rise to the observed relation over about two orders of magnitude 
in Up, when observed off-axis along different sight lines, provided the opening angle and 
angular width of the jet are larger than 1/T, where T is the Lorentz factor of the emitting 
material (see appendix for an approximate analytic derivation). The effect is not necessarily 
specific to annular jet; it applies to any geometry in which the emission along the line of sight 
is comprised of contribution from an extended source (e.g., Toma, et al. 2005). However, 
there are a priori reasons to favor an annular geometry (see below). EL04 also calculated 
the rate distribution of observed peak energies u p , and have shown that the relative number 
of XRFs observed is compatible with the model, provided the angular width of the emitting 
region is not too large. (Off-axis effects have been considered also by Yamazaki et al. [2004, 
and references therein]. However, they assumed an ad hoc intrinsic relation between E iso 
and v v designed to fit the observations.) If off-axis effects are affecting v v and Ei so then the 
various quantities one derive from GRB observations need to have this effect folded in. For 
example, Eichler & Jontof-Hutter (2005) have recalculated the average prompt gamma ray 
to blast wave energy ratio and found it to be considerably higher than previous estimates, 
simply because the prompt gamma rays are more affected by viewing angle effects than the 
blast wave energy as inferred from the afterglow. The normalization of the Ei SO — v v curve 
produced by the model of EL04 depends on the amount of energy per solid angle ejected by 
the source, which for a standard energy release is inversely proportional to the solid angle 
subtended by the emitting region. Evidently, any spread in this parameter in a sample of 
sources should be reflected as a scatter in the resultant correlation. Under the assumption of 
a standard energy output, this scatter can be reduced essentially by correcting for beaming. 
If an absolute way can be found to correct for beaming, then the model naively predicts the 
collimation-corrected energy to satisfy the same relation as the isotropic equivalent energy, 
but with a different normalization and virtually no scatter. 

An attempt to derive a relation between the collimation-corrected GRB energy, denoted 
here by E ltapp for future purposes, and the peak energy has been made recently by GGL04. 
These authors considered a sample of 24 GRBs with measured redshift and peak energy, for 
which the break time of the afterglow light curve, tt, rea k, is well constrained. They found 
indeed a tight relation between E 1APP and the observed peak energy z/ p , but with a different 
slope: E ljapp cx v p ' 7 . At first sight this seems to be in conflict with the prediction of the model 
discussed above. However, careful examination reveals that this result is in fact consistent 
with the hypothesis that the Amati relation is due to viewing angle effects, and that the 
change in slope is primarily due to a bias in the beaming correction used by GGL04. 

In order to correct for beaming GGL04 estimate the opening angle of the emitting jet 
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for each source in their sample, using the two observables, E iso and threat, and the relation: 

/ , n 3/8 / \ 1/8 

e(t break ,E lso ) = 0A6l( t -^) . (2) 

td is the break time measured in days, ?? 7 is the radiative efficiency, and n the 
density of the circumburst ambient medium, assumed to be constant. Since the radiative 
efficiency is unknown, they invoked a single value for all sources of i] 7 = 0.2. For the 
circumburst density they adopted the value n = 3 cm -3 for all sources, except five for which 
estimates for n are available in the literature. The collimation-corrected energy was then 
taken to be: E ltapp = (1 — cos 6)E iso ~ (9 2 /2)E iso . Since the opening angle given by eq. (2) 
depends on E iso , it is clear that it may be subject to biases in cases where E iso is influenced 
by effects other than the spread in the opening angle. This is particularly true for sources 
observed off-axis, for which eq. (2) underestimates the beaming factor. Furthermore, the 
break time of the afterglow emission may appear longer for off-axis observers (Granot et al. 
2002), leading to an even larger bias. The latter effect is typically small, except for the very 
soft sources (i.e., those observed at the largest viewing angles). In the following we provide 
a quantitative treatment of these biases. 

The model outlined in EL04 assumes a uniform, axisymmetric jet of opening angle 9 2 , 
with intrinsic spectral peak at v*, and with a hole of angular size Q\ cut out of it. This 
symmetry was chosen strictly for convenience, and it can be easily seen that the scaling de- 
rived in EL04 applies to more complicated geometries. The peak frequency and the observed 
isotropic energy were calculated numerically for different viewing angles inside the hole and 
outside the jet, and it was found that a relation Ei SO oc v" holds over a wide range of u p , 
with a laying in the range between 2 and 3, depending on source parameters. In particular, 
a ~ 2 when 9 2 and A9 = 9 2 — Q\ are larger than T -1 . The analysis of EL04 does not account 
for a possible spread in source parameters. In particular, the normalization of the E iso — u p 
relation produced by the model, depends on the angular distribution of the released energy. 
To illustrate how this might affect the observed E iso — v p relation, we consider a sample of 
sources having a universal spectrum and a standard energy output, E 7 , but a range of open- 
ing angles. Denote by A b = 47r/Af2 the corresponding beaming factor, where AQ is the solid 
angle subtends by the gamma-ray emission region. For a symmetric (double sided), hollow jet 
we have AQ = 4vr f£ sin 9d9 = 4n(cos9 1 -cos# 2 ), and A b = (cos0i - cos #2)^ ^ 2/(91-61). 
The isotropic equivalent energy measured for a source observed at some angle corresponding 
to an observed peak energy u p , can then be expressed as, 

Eiso = A b E^{y p j v p ^ max ) . (3) 

Here v P)Tnax = 2Yu* denotes the peak frequency of a source observed head-on. Consequently, 
a sample of sources with fixed E 1 and a range of beaming factors A b would form a family 
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of parallel lines in the E iso — v v plane, with each line corresponding to a subset of sources 
with beaming factors in the interval + dA^). It is tempting to interpret figure 2 in 

Ghirlanda et al. (2005) and figure 5 in Bosnjak et al. (2005) as such. The sample studied 
originally by Amati (2002), that contained particularly bright bursts, may consist mainly of 
the subset of sources with the largest beaming factors, and therefore defines a boundary in 
the E iso — v v plane, as suggested by Nakar & Piran (2004). 

There are several reasons why the GRB fireball might be shaped like a thick annu- 
lus. First, the collimation of supersonic material can itself produce an annular shaped, high 
entropy jet (Eichler 1982, Levinson & Eichler 2000; Begelman & Blandford, private com- 
munication) because the jet material tends to accumulate into a shock-compressed layer at 
the confining walls. Second, the soft gamma ray emission is likely to originate from that 
part of the jet that is baryon loaded, and the baryon loading may come from the periphery. 
In the particular case of baryon loading by neutron leakage from the walls of a confining 
wind or stellar envelope (Levinson & Eichler 2003), the neutrons are quickly charged by 
collisions near the walls, before they can penetrate to the center. It was shown in the above 
reference that the annular region that is significantly loaded can have a solid angle that is 
not too much less than the inner hollow region. Third, gamma-rays emitted by a compact 
photosphere (Eichler & Levinson, 2000) can mostly impinge on the baryon rich periphery of 
the outflow and be dragged along by it, and thus be concentrated into a an annulus with 
a smaller total solid angle than that into which they were originally emitted. The inner 
core may be comprised of Poynting flux that contribute very little to the soft gamma-ray 
emission, but nonetheless carry a considerable fraction of the ejected energy. Thus, even 
though the gamma-ray emission region may be annular, the total energy carried by the out- 
flow is likely distributed more uniformly inside the cone. Now, the structure of the blast 
wave driven by the collision of the fireball and the surrounding gas depends mainly on the 
angular distribution of total energy of the piston. It is therefore reasonable to assume that 
at small viewing angles the break in the afterglow light curve is associated with the opening 
angle of the outer cone. At large viewing angles, the break time of the afterglow emission 
may be altered. The exact shape of the afterglow light curve as viewed by off-axis observers 
may depend on details. For sight lines outside the jet, one might naively expect a break 
when the Lorentz factor drops to T ~ i/0 n , where 9 n is the viewing angle measured with 
respect to the jet axis (Granot et al., 2002). The angular separation corresponding to an 
observed peak energy v v is given by: 9 n — 9 2 = ~£~ x ^v v ,maxlv p — 2) 1 / 2 . The distribution 
of opening angles inferred from the achromatic breaks of the afterglow emission peaks at 
6 ~ 0.1 (Frail et al. 2001; Guetta et al. 2005; Ghirlanda et al. 2005). Adopting this value 
for 9 2 we find that 9 n — 9 2 < 9 2 if v v , ma xlv v < T 2 /200. For reasonable values of the Lorentz 
factor T this condition is satisfied essentially for all sources in the GGL04 sample, and we 
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therefore anticipate the break time not to be significantly altered by viewing angle effects. 
For observers looking down the hole the break time may appear unaltered. Let us denote 
by q( v p) = tbreaki^pj/tbreaki^max) > 1 the ratio of break times measured by off-axis and 
head-on observers. Suppose that for a source observed head-on, the opening angle obtained 
by substituting the observed energy, Ei SO = A b E 1 , and break time, tbreak, into eq. (2) is 
roughly 9 2 , viz., 8(t break { ^p,max)i A-bE^) — 9 2 . Then for the same source observed off axis the 
inferred opening angle would be, 9(t break , E iso ) ~ 9 2 q 3/8 (v p /v Ptmax )- a/8 , where eqs (2) and 
(3) have been used. The apparent collimation-corrected GRB energy, E^^ app , is then given 
by 

E^^app T\ ^-bE^iyp j ' Vp^max) — Z T7, , Q \y Ey(Up/Up tmax ) . (4) 
I 1 — (01/ 02 J 

For a = 2 and q — 1, the last equation yields: 

E^^app oc z-'p' 5 , (5) 

in agreement with the result obtained by GGL04, -E 7j(W cx z/i- 416±0 - 09 . We expect the index 
to be somewhat smaller than 1.5, owing to break time effects contributed by the softest 
sources, those with v p « u p , max , for which the ratio q may already be large enough to affect 
the correlation. As seen from eq. (4) the value of E 1APP predicted by the model depends in 
addition on the ratio 9\j9 2 . This parameter is uncertain in the present picture. A spread 
in this parameter will contribute a scatter in the observed E lyapp — v p relation. The inverse 
correlation between E iso and jet opening angle (Frail et al. 2001; van Putten & Regimbau 
2003), seems to imply that the annulus is reasonably thick, such that it subtends a significant 
fraction of the solid angle subtended by its outer cone. In this case the scatter is expected 
to be rather small. For the range of parameters adopted in EL04 the values of the coefficient 
1/[1 — (9i/9 2 ) 2 ] in eq. (4) vary between 1 and 1.5. 

The above theoretical explanation is consistent with the observed correlation between v p 
and the inferred opening angle (Eichler & Jontof-Hutter, 2005), where the inferred opening 
angle shows the tendency to increase with decreasing v p . The point is that the inferred 
opening angle is weakly biased towards larger values when E iso is underestimated due to off 
angle viewing (see eq. 2). 

Whether the observed relations discussed above hold for the entire population of GRBs 
or are the result of some selection effects is at present under debate. Several groups have 
attempted to perform consistency checks for large samples of BATSE sources with known 
fluence and well determined peak energy. Nakar & Piran (2004) analyzed trajectories in 
the E iso — v p plane obtained for sources with a measured fluence and observed peak energy, 
v p ,obs = v p /{l + z), by varying the redshift. They then define a source to be an outlier if 
the minimum distance between its trajectory and the curve representing the Amati relation 
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exceeds a certain value. By applying this test to a sample of bright bursts from (Band et 
al. 1993; Jimenez et al. 2000) they concluded that about 50% of the sources in their sample 
are outliers. However, they also concluded that the line representing the Amati relation 
defines a boundary of the region in the E iso — v v plane which is populated with sources. 
In other words, all outliers should be dim, hard bursts. As argued above, this boundary 
may represent the subsample of GRBs with the smallest opening angles. Band & Preece 
(2005) extended this work to a larger sample of BATSE sources and tested in addition the 
Ghirlanda correlation. They find an even larger number of sources to be inconsistent with the 
Amati relation. The number of outliers to the Ghirlanda relation depends on the assumed 
distribution of beaming factors, and appears to be much smaller, at least for certain choices 
of beaming correction. This better agreement may be due to the smaller scatter in the 
Ghirlanda correlation, although we are aware of the large uncertainty in the distribution of 
opening angles that might affect the significance of this result. Different consistency checks 
have been performed subsequently by two other groups. Ghirlanda et al. (2005) analyzed 
a sample of 442 bursts with what they term "pseudo redshifts" which have been estimated 
using the lag-luminosity relation (Norris et al. 2000; Norris 2002). Bosnjak et al. (2005) 
tested the consistency of the Amati relation with the fluence distribution of bright BATSE 
bursts assuming that the GRB population follows the star formation rate. Both groups 
concluded that their samples are consistent with the Ei SO — v p relation discovered by Amati, 
but with a larger scatter than originally found (but cf., Nakar & Piran, 2005). 

Here we are arguing that the scatter is to be expected. Clearly, if there are two effects 
that are creating two separate correlations, a study of just one of them will find scatter. 
As suggested above, both opening angle effects (Frail et al. 2001) and viewing angle effects 
(EL04) are each creating their own correlation: the opening angle effect would create scatter 
in E iso even for head-on viewers and even if there were no scatter in i? 7 , and a viewing 
angle effect would create scatter in the observed E iso even if there were no scatter in E 1 
or in E^/9 2 . In addition, extreme outliers such as GRB 980425 can be attributed to a 
relatively large angle scattering off material with a relatively modest Lorentz factor (Eichler 
& Levinson 1999; Nakamura, 1998). Assuming that this component is weak it can probably 
only be seen from very nearby GRBs. 

In conclusion, a population of beamed sources with a standard energy output and a 
universal intrinsic spectrum, viewed along different sight lines, can explain the Amati and 
Ghirlanda relations. The small scatter in the peak energy - collimation corrected energy 
correlation reported by GGL04 indicates that the scatter in the E iso — v v relation is predom- 
inantly due to a spread in the solid angle subtended by the gamma-ray emission region, and 
that the total energy output in gamma-rays of long GRBs is narrowly clustered. The total 
energy released as gamma rays should equal roughly the corrected energy E ltapp of the hard- 
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est source in the GGL04 sample, E 1 ~ 2 x 10 51 ergs. Consistency checks of the Amati and 
Ghirlanda relations on large samples should allow for the potentially large scatter produced 
by the opening angle effect and for biases in the beaming correction factor. Larger samples 
of sources with measured redshifts, particularly dim soft bursts, are ultimately needed to 
confirm those relations. The original argument of Frail et al. (2001), that the scatter in 
Ei SO is greatly reduced when the beaming correction they make is invoked may have naively 
appeared to have been confounded by the large spread in apparent E iso that remained. Here, 
however, we have noted that the observable difference between the slopes of the Ghirlanda 
and Amati relations in fact provides an interesting confirmation to the afterglow theory in- 
voked by Frail et al. (2001). In particular, the noticable differences between the Amati and 
Ghirlanda slopes confirms the very weak dependence of opening angle on Ei so , which may 
have otherwise gone without direct observational confirmation. 

This research was supported by the Israel-US Binational Science Foundation, an Israel 
Science Foundation Center of Excellence Award, and the Arnow Chair of Theoretical Physics. 



A. Analytic derivation of the E iso — v v relation 

Consider the emission from an annular jet centered around the 2-axis, and having an 
opening angle 9 and angular width A9. Denote by (5 the directions of emitting fluid elements 
(assumed to have a fixed Lorentz factor Y) and by h the sight line direction, respectively, and 
define cos 9p — f3 • z, and cos 9 n = n- z. Then h • f3 — cos 9p cos 9 n + sin 9p sin 9 n cos 0, where <fi 
is the azimuthal angle. Assuming that each fluid element emits isotropically in its rest frame 
a total energy E', then its contribution to the overall energy emitted into a solid angle dQ n 
around the sight line direction n, as measured in the Lab frame, is: dE(n,/3) = (E'V 3 )dQ n , 
where D = 1/T(1 — f3 ■ n) is the corresponding Doppler factor. The fluence along the line of 
sight is given by the integral 



E = 



[ dE ^^ dQp oc [ P 3 sinVV0- (Al) 



Define A = 9 n — 9/3, and assume small angle approximation, viz., A# < #o << 1 and A < #0i 
the Doppler factor can be expressed as: 

1 2 
P ^ A ' X ^ = r[l -/?cosA + /5sin^sin^(l -cos0)] ~ r[r- 2 + A 2 (l + X 2 )} ' ^ 

where x 2 — 2 sin 9p sin 9 n (l — cos0)/A 2 . Define rf = 2 sin 9p sin 9 n /A 2 , and 

F(A) = f o + xw-xv^ - (A3) 



- 8- 



equation (Al) can be written, for angular separations A 2 » T 



as 




8 /sin(0 n -A) 



sin 6, 



F(A)dA, 



(A4) 



with A min = 9 n — 6 , A max = A min + A9. In the limit considered here, viz., A << O < 9 n , 
we have i] >> 1, so that F(A) ~ / °°(1 + x 2 )~ 3 dx = 37r/16, and sin(#„ - A)/sin# n = 
1 — O(A/0 n ). Equation (A4) reduces to Ei SO oc r~ 3 A~| n . For a jet emitting an intrinsic 
single universal spectrum with a peak frequency u*, the observed spectral peak will be located 
at Up = u*V(A min , x = 0) = 2u*/ (rA^ in ) for sight lines outside the jet (EL04), whereby we 
obtain E iso oc r~ 1 (z/ p /z/*) 2 . Similar result can be readily derived for sight lines inside the 
hole. Numerical integration of eq. (Al) for a wide range of parameters and viewing angles 
is presented in EL04. 
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